Stabilized vortex solitons in layered Kerr media 



Caspar D. Montesinos and Victor M. Perez-Garcfa 
Departamento de Matemdticas, Escuela Tecnica Superior de Ingenieros Industriales, 
Universidad de Castilla-La Mancha, 13071 Ciudad Real, Spain 



: 

o . 
o ■ 

(N ■ 



(N 

in 

PL, 

d 



(N 
> 

in 
o 
in 
o 

o 



> 

x 



Humberto Michinel and Jose R. Salgueiro 

Area de Optica, Facultade de Ciencias de Ourense, 
Umversidade de Vigo, As Lagoas s/n, Ourense, ES-32005 Spain. 
(Dated: February 8, 2008) 

In this letter we demonstrate the possibility of stabilizing beams with angular momentum propa- 
gating in Kerr media. Large propagation distances without filamentation can be achieved in layered 
media with alternating focusing and defocusing nonlinearities. Stronger stabilization can be ob- 
tained with the addition of an incoherent beam. 
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Introduction.- Vortices have been a source of fasci- 
nation since the works of Empedocles, Aristotle and 
Descartes, who tried to explain the formation of the 
Earth, its gravity and the dynamics of the solar sys- 
tem as due to primordial cosmic vortices. Many inter- 
esting problems related to vortices are open in different 
fields such as: fluid mechanics, high T c superconductivity, 
superfluidity, light propagation, Bose-Einstein condensa- 
tion (BEC), cosmology, biosciences, or solid state physics 

Quia. 

In wave mechanics a vortex is a screw phase disloca- 
tion, or defect where the amplitude of the field van- 
ishes. The phase around the singularity has an integer 
number of windings, £, which plays the role of an angular 
momentum. For fields with non- vanishing boundary con- 
ditions, this number is a conserved quantity and governs 
the interactions between vortices as if they were endowed 
with electrostatic charges. Thus, I is usually called the 
"topological charge" of the defect. 

In Optics there has been a strong interest on the so 
called "vortex solitons" , i. e. robust distributions of light 
of vortex type in which nonlinearity could compensate 
diffraction leading to stationary propagation. However, 
in self-focusing Kerr media, a finite size beam contain- 
ing a vortex always destabilizes and forms a filamentary 
structure |(|. This also stands for saturable self- focusing 
nonlinearities Q. Vortex solitons have been studied in 
many other different optical systems (see e.g. the review 
H) and in most realistic cases they tend to be unstable. 

In this paper we propose to use layered Kerr media, 
which are self-focusing in average, to obtain stable prop- 
agation of vortex solitons up to very long distances. Our 
ideas are also extended to the field of matter waves. 

Stabilized solitons.- In Ref. [2j an idea to prevent 
collapse was proposed based on the modulation of the 
Kerr coefficient of an optical material along the propaga- 
tion direction. In that way, a propagating beam focuses 
and expands in alternating regions and becomes stabi- 
lized in average. This idea has been further explored in 
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The propagation of a paraxial monochromatic beam in 
a Kerr medium is modeled by equations of the type (in 
adimcnsional units) 
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where u(x,y,z) : R 2 x M. + — > C is the slowly varying am- 
plitude of the beam envelope, A = d 2 /dx 2 + d 2 /dy 2 , and 
g(z) is a periodic function accounting for the modulation 
of the nonlinearity. It is well known that, if gf is constant, 
there is a stationary radially symmetric solution of Eq. 
JU (the so-called Townes soliton): u(r, z) — <&o{r)e lXz . 
This solution is unstable since any generic slight pertur- 
bation of the initial condition will yield either collapse 
or spread of the distribution. In previous works, it has 
been shown that the structure which arises when the non- 
linearity is modulated as described above is a stabilized 
Townes soliton (STS) 00. 

Partially stabilized vortex solitons.- However, Townes 
solitons are not the only stationary solutions of Eq. (JTJ 
for constant g. There also exist vortex- type solutions of 
the form u(r,z) — ^g(r)e iee e lXz which are unstable as 
well. Therefore, one could naively expect that the same 
stabilization mechanism proposed in Ref. could be 
applied to achieve stabilization of these solutions, i.e. to 
induce alternative expanding and squeezing cycles of the 
vortex by a periodic modulation of the nonlinear coeffi- 
cient. To study further this possibility we have computed 
numerically the profiles of vortex solutions of Eq. (|T|) 
corresponding to a constant value of g v = —24.15 (the 
critical one for vortex solutions), by using a standard 
shooting method. We have then studied the evolution 
of this stationary solution numerically using a special 
pseudospectral method described in Ref. [1J| in differ- 
ent situations. Obviously if we compute the evolution of 
this vortex solution for constant g > g v we would have 
expansion while for g = go < g v it would collapse. In 
Fig. nja,b,c) we plot the collapse of a vortex with su- 
percritical value of go = — 8ir, while in Fig. ^d,e,f) we 
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FIG. 1: [Color online] Evolution of initial data of vortex type 
obtained as a stationary solution of Eq. Q for g v — —24.15 
under different parameter variations, (a-c) Evolution for con- 
stant g — —8n. (d-f) Evolution with modulated nonlinearity 
g(z) = — 8n + 207rcos40,z. (g) Isosurface plot of u(x,y,z) of 
the same simulation as in (d-f) showing the development of 
the instability, (h-j) Evolution with modulated nonlinearity 
g{z) = -27 + 20ttcos40z. 



present some snapshots of the evolution of the vortex af- 
ter the addition of a stabilizing term to the nonlinearity 
g(z) = go+gi cos Qz for g = —8n, g\ = 20n 7 il = 40. We 
can see how the periodic modulation of the nonlinearity 
retards the filamentation of the vortex, which now prop- 
agates for a longer distance before breaking into several 
stabilized solitons. The splitting length and number of 
outgoing stabilized solitons depend on the specific param- 
eters of the nonlinearity. For instance choosing a different 
set of parameters go = —27, g\ = 20tt, O = 40 increases 
the stability of the vortex and decreases the number of 
filaments. Since each emerging beam is close to an STS 
the excess energy is eliminated in the form of radiation 
which is removed by the absorbing boundary conditions 
of our numerical scheme. In this paper we have chosen a 
smooth form for the modulation g(z) but similar results 
are obtained when g(z) is taking as piecewise constant. 

We have made an extensive search in the parameter 
space for modulations of the form g(z) — go + gi cosfiz 
and we have not been able to find any parameter combi- 



nation allowing indefinite stabilization of the vortex soli- 
ton. Concerning finite-dimensional reduced models for 
the evolution of the effective width of the solutions, such 
as those successfully used for nodeless beams in Refs. 
[ToL ITU Jl^L Hii ]. we must stress that as pointed out in 
Ref. [l0( these formulations do not reflect correctly the 
dynamics and instabilities of vortex solutions. 

Vector systems.- From the previous analysis it seems 
that a vortex can only be partially stabilized in the frame- 
work of Eq. QJ, i.e. in scalar systems. Recent works 
point out the important fact that the incoherent interac- 
tion of two components could provide, in saturable me- 
dia, an effective wave guid e for the vortex, leading to a 
more stable behavior [Ifj. Following this idea we con- 
sider now a vector two-component system with Kerr in- 
teractions of the form: 

du 1 

i-^- = --Ami + g(z) (an|ui| 2 + ai2|u 2 | 2 ) wi,(2a) 



dz 
.du 2 

dz 



Au 2 + g(z) (a 2 i|ui| 2 + a 2 2\u 2 \ 2 ) u 2 ,(2b) 



where ajk G K are the nonlinear coupling coefficients and 
g(z) accounts for the modulation of the nonlinearity. We 
will denote Ij — L 2 \iij\ dxdy. Although this system is 
conservative, in the numerical simulations to be shown 
later we incorporate absorbing boundary conditions in 
order to get rid of the radiation. Therefore, in practice, 
there will be a decrease of Ij during the propagation. 

Eqs. J2J are an extension of the Manakov system [J^ 
to two transverse dimensions. Among other situations 
these equations model the propagation of two circularly 
polarized beams with opposite polarizations leading to 
specific factors an = a 22 = 1, a\ 2 = a 2 \ = 2. In the con- 
text of BEC, these equations (with an additional trapping 
term) describe the dynamics of multicomponent quasi- 
two dimensional condensates, uj being the wavefunctions 
for the atomic species involved. The formation of vec- 
tor solitons composed of appropriate fractions of Townes 
states has been studied in Ref. [l^ . 

Our idea is to choose g(z) to get the Townes soliton 
Ui stabilized. As the coupling terms in Eq. J2J) would 
provide an effective waveguide for u 2 , it seems reasonable 
that, when I2 <C h, the guiding effect will dominate over 
self-interaction and the vortex could become stabilized. 

Limit of small u 2 .- Let us first consider the case of 
constant g and I 2 <C I\ so that Eqs. J5J) become: 



.du\ 1 . l2 

1—— ~ —-Aui+ga n \ui\ m, 
oz 2 
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(3a) 
(3b) 



Taking Ui(r, z) = & Q (r)e lXoz , then Eq. if3b)l is a linear 
two-dimensional Schrodinger problem in which the role 
of the potential is played by |$o| 2 - Following Ref. 0] 
we can bound the number of ^-wave bound states in this 
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FIG. 2: [Color online] Solutions of Eqs. |(2J for initial data 
ui(r, 0) = $o('"),'"2 = a^(r)e l9 with a = 0.1 in a grid of 
810 x 810 points on the spatial region [-40,40] showing stable 
(but dissipative due to the effect of radiation) propagation of 
the vortex in the full range z £ [0, 500]. (a) Isosurface plot of 
U2(x,y,z) spanning all the propagation range, (b) Evolution 
of the norms I^ 2 (z) , I^ 2 (z) and of the amplitudes Ai(z) — 
max( IiS ) |ui|, A 2 (z) — max^^j |ii2| of both components. 

potential, by Ng < (ai 2 g/£) J R+ r\ui\ 2 dr . Vortex-type 
solutions with smallest topological charge are those with 
£ = 1. Thus, taking into account that for a Townes 
soliton gl\ — 0.931 we get Nz—\ < 0.931ai2- Therefore, 
we can expect the existence of a unique £ = 1 stationary 
vortex solution of Eq. (|3b|) in the case 0,21 = 2. This 
will be denoted hereafter as it 2 — V(r)e te e lX ^ z . Let us 
notice that for £ > 2 we get always Ng>2 < 1 thus ruling 
out the possibility of obtaining higher order vortices. We 
have numerically found the profile of this vortex solution 
V(r) by using a standard shooting method. 

Choosing an appropriate modulation for g allows us to 
stabilize u\. For Eq. I)3b|l . the potential -the stabilized 
Townes soliton- oscillates with a fast frequency of the 
order of f2 and another slower one of dynamical origin 
|l3|. In this linear case we can apply to it 2 the quantum 
mechanical theory of fast perturbations |l9| to account 
for the effect of the fast modulation on the vortex. The 
main result of this theory is that the vortex will remain 
unaffected by the fast perturbation in the potential pro- 
vided the modulation period T satisfies A U2 H <C 1/T. 
In our case H — ^ J Q T [— + g(t)a,2i \ui | 2 ] and this 
inequality imposes T <C 1.12 which requires f2 3> 5.6. 
The STS oscillation of lower frequency will induce the 
same modulation in the vortex. To verify our previ- 
ous considerations we have simulated Eqs. Q with ini- 
tial conditions ui(r, 0) = &o{r) 7 U2(r, 0, 0) = V{r)e l6 and 
g{z) = — 27r+87rcos40z. Full stabilization of the vortex is 
observed, its oscillations following the pattern predicted 
above, i.e. there is only a residual fast oscillation in the 
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FIG. 3: [Color online] Same as Fig. Hfor a = 0.32. (a)-(e) 
Pseudocolor plots of u 2 {x,y,z) for z = 0,200,295,305,450. 
(f) Isosurface plot of U2(x,y,z) spanning all the propagation 
range, (g) Details of the region in which the vortex spirals 
out of U2 for z £ [280,319]. (h) Evolution of the norms of 
both components I^ 2 (z), I^ 2 (z) showing the readjustment 
of the norms after the vortex is ejected and a stabilized vector 
soliton is formed, (i) Amplitudes Ai(z) = max( IiS ) |ui| and 
A 2 {z) = max( 1|S ) |«2| of both components. 



vortex component and its slow oscillation follows that of 
the stabilized Townes soliton in u\. 

Fully nonlinear regime.- To account for the case of fi- 
nite 112 we have studied numerically the solutions of Eqs. 
@ taking as initial data ui(r, 0) = $o( r ), U2(r, 8, 0) = 
aV(r)e and g(z) = —2tt + 87rcos40z. We have re- 
peated most simulations starting with the full stationary 
solutions of Eqs. @ and found similar results. 

For small values of a (e.g. a — 0.1 as in Fig. [2J the 
vortex is fully stabilized up to the maximum propagation 
distances studied. However, a continuous loss of energy is 
observed during propagation. We think that this power 
damping is due to radiation emitted by the vortex and it 
is related to the continuous background oscillations of the 
stabilizing Townes beam. We have tested these numeri- 
cal simulations both by the spectral numerical scheme of 
Ref. ^4|], for which radiation is eliminated by means of 
an absorbing potential and by a finite-difference Crank- 
Nicholson type scheme with transparent boundary con- 
ditions with the same results. 
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For larger values of a (e.g. a = 0.32 as in Fig. |3J), the 
vortex destabilizes at long propagation distances due to 
the effect of nonlinear interactions between the guiding 
Townes soliton and the vortex. We can see that, although 
the perturbation is not small, the vortex propagates for 
very long distances of about 300 propagation units (com- 
pare this with the results shown in Fig. In the region 
of stable propagation the vortex amplitude slowly decays 
due to the emission of radiation until a stabilized vector 
soliton is formed. In this process, both components emit 
radiation to readjust their norms to satisfy the relation 
h+h = I T o W nes M From Fig. Hi) it is clear that 
the oscillations of the vortex amplitude basically contain 
only the slower frequency and that our previous argu- 
ments apply here. The addition of noise to the initial 
data of 1% in amplitude triggers the instability faster but 
even in that case the vortex propagates for more than 125 
adimensional units before the instability sets in. 

Finally let us briefly comment that for a = 0(1) there 
is a different branch of bound states of Eqs. (JSJ) in the 
fully nonlinear regime, corresponding to thin Townes soli- 
tons localized near the bottom of the vortex solutions. 
This soliton act as a pinning potential for the vortex but 
it is not enough to stabilize vortex- type solutions. 

Applications to matter waves.- The previous results 
have also implications in the field of matter waves be- 
cause of the close analogy of Eqs. J2J) with the equations 
of evolution of a multi-component Bose-Einstein conden- 
sate in the mean field approximation. The analysis of 
vortices in dilute-gas BECs has been a very hot topic in 
the last years, specially after their experimental genera- 
tion with different setups Q. In multicomponent BEC 
systems the interaction coefficients are proportional 
to the respective scattering lengths and the effective two- 
dimensionality can be achieved by confining the conden- 
sate tightly along one specific direction. The condition 
Ni < 0.931ai2 has relevance since it imposes restrictions 
to the atomic species which can be used to trap a vortex. 
For instance the cross- interaction coefficient in multicom- 
ponent condensates made of different hyperfine levels of 
87 Rb does not satisfy this condition. However, bosonic 
K-Rb mixtures such as the one described in Refs. [2(1 EH ] 
could be used because of the large scattering length of 
the collisions K-Rb. In this scenario the results pre- 
sented here could be extended to tightly confined BEC 
systems. Our predictions would imply the existence of 
self-supported vortex solitons which could be generated 
using Feschbach resonance management techniques. 

In conclusion, we have shown stabilization of vortex 
solitons for long propagation distances in stratified Kerr 



media by control of the nonlinear coefficient and stronger 
stabilization by the use of two combined beams. Our 
results can also be used to stabilize vortices in certain 
types of multicomponent Bose-Einstein condensates. 
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